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Abstract 

This article attempts to delineate the roles played by non-dynamical background 
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tensors generated from a diffeomorphism invariant action density. An intrinsic coor- 
dinate independent approach puts into perspective a number of spurious arguments 
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Hilbert stress-energy-momentum tensor derived in the context of Einstein's theory 
of general relativity. Emphasis is placed on the role played by non-dynamical back- 
ground (phenomenological) structures that discriminate between properties of these 
tensors particularly in the context of electrodynamics in media. These tensors are 
used to construct conservation laws in the presence of Killing Lie-symmetric back- 
ground fields. 
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1 Introduction 



The venerable aim of the variational calculus as a tool for deriving from a single princi- 
ple the observed laws of physics remains an attractive one. It synthesises conservation 
laws from local symmetries and offers a route for finding unified schemes that under- 
pin the experimental sciences. However the epistemology that arises in this approach 
has, on occasion, lead to unnecessary confusion when employed in the discussion of sys- 
tems that require phenomenological input. In particular this remains true in the context 
of electrodynamics in macroscopic material where the indiscriminate use of the stress- 
energy-momentum concept in classical physics has lead to different notions of quantum 
electrodynamics in such media. The calculus often enables observed local and global sym- 
metries of nature to be accommodated in a natural manner and can be used to generate 
conserved currents associated with local symmetries. Furthermore the local dynamics 
of coupled systems follows from local extrema of a functional constructed from a single 
diffeomorphism invariant action integral S = J M A T , for some diffeomorphism invariant 
differential 4- form A T on spacetime Ai. If one chooses a nowhere vanishing 4- form Q 
on spacetime and writes A T = SQ then S is called a scalar action density (relative to 
Q). For any covariant metric g on Ai a natural choice is Vt = ^ ] /\detg\ d 4 x = *1 where 
* is the Hodge map associated with g. An alternative choice, given a coordinate system 
(ar, . . . , x 4 ) is the tensor density Q = dx l A • • • A dx . Such a choice implies that S then 
transforms like | det t7 1 x / 2 under a local change of coordinates^. 

Traditional symmetry analysis based on the pioneering work of Noether exploits such 
densities often in a manifestly coordinate dependent manner. Modern approaches exploit 
an intrinsic jet-bundle formulation. Between these extremes are ad-hoc formulations that 
have led to intense debates about the role of alternative stress-energy-momentum tensors 
used to describe the interchange of energy and momentum between interacting systems in 
continuous macroscopic media. In particular, conserved Noether currents have historically 
been constructed from 3-forms on space and their associated Lagrangian scalar densities. 
Notions of stress and power in continuum mechanics became incorporated into the Lorentz 
covariant formulation of special relativity. With the advent of Einstein's general theory 
of relativity the spacetime symmetric metric tensor field becomes a dynamical variable 
and an alternative stress-energy-momentum tensor can be defined through the variational 
derivative of the action 4-form on spacetime (for gravity and matter) with respect to this 
metric tensorJl The issue that then arises is how best to identify conserved quantities 
constructed from different models that can be put into a variational formulation and 
relate them to different choices of stress-energy-momentum tensor. 

Conserved quantities can be generated from closed 3-forms in spacetime. These often 
arise as a consequence of some symmetry of an element or elements contributing to the 
structure of the action for the model. However even in the absence of such symmetries 
there are powerful relations that arise from applying the variational approach to action 
functionals that are (locally) invariant under transformations of the fields induced by 
arbitrary (local) spacetime diffeomorphisms. Such functionals are readily constructed in 

J If U = d/dx 1 generates time translations then ijjh T may be called a Lagrangian 3-form [3] relative 
to U. 

2 Matter in this sense includes electromagnetism as well as other matter fields. Modified theories of 
gravity, e.g. with scalar fields coupled directly to spacetime curvature, offer no unique identification of a 
"matter" action. 
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terms of coordinate invariants made by contracting tensor fields of various degrees or 
tensor "densities" of different weights. It is not necessary that all such quantities in 
this construction be dynamical, i.e. subject to variational field equations. However it is 
in these circumstances that there arise differences when one compares the consequences 
deduced from different choices for the description of stress, energy, and momentum in the 
presence of such quantities. 

The (covariant) variational approach is often restricted to closed (non-dissipative) 
systems where dynamical equations for the dynamical field variables arise by finding local 
extrema of a total action S under their variations, where A T is a 4-form on a spacetime 
manifold. Such forms belong to a class, members of which describe the same classical 
physics. For example members that differ by an exact 4-form with compact support yield 
the same variational equations. 

Since the description of gravitation is given in terms of a geometry of spacetime, any 
collection of dynamical field variables {Cl_ = {&> • • • > Cq} that are n °f P ar t °f the geometry 
may be assigned the status of "matter" o In Einstein's theory, gravitation arises from a 
pseudo-Riemannian spacetime geometry based entirely on a dynamical spacetime metric 
tensor field g. In the absence of matter, the dynamics of g is given by the Einstein-Hilbert 
action J M A Em . The dynamics of (non-spinorial) matter minimally coupled to gravitation 
is described by a matter action 4-form A m that is independent of derivatives of g. 

In order to eliminate use of the jet-bundle language and simplify the analysis, in the 
following the matter action 4-form A m will be restricted to depend on g, {Z}, {£} and 
{d£}, where a collection of prescribed non-dynamical background tensor fields {Z} = 
{Zi, . . . , Zn} has been included. Such fields play no variational role in fixing the local 
extrema of S but (if present) play a crucial role in determining the consequences that 
follow from diffeomorphism invariance. Then the total action 4-form is A T = A Em + A m . 

It is sometimes useful to define matter subsystems that are described by sub-actions 
f M A s where A m = X] s ^ s - Clearly this notion of a subsystem is defined relative to a 
particular decomposition of A m and extrema of S will not in general coincide with extrema 
of J M (A Em + A s ). When such extrema do coincide one may argue that the subsystem 
described by A s decouples from the system described by A T , i.e. it becomes a closed 
subsystem. This is rarely the case in systems interacting with dynamical gravitation. 
Subsystems that are not closed are termed open. 

Einstein's gravitational field equations in the presence of matter lead one to iden- 
tify the variational derivative of A m with respect to g (see below) with the symmetric 
"stress-energy-momentum" tensor T H (g, {Z}, {£}, {<i£}) associated with matter. This 
terminology is natural, given its historic connection with Newtonian concepts but ar- 
guably misleading in a broader context where symmetries associated with space and time 
translation are absent. The tensor T H relies for its definition on A m depending on the 
spacetime metric. On manifolds where gravitation is irrelevant (where g is regarded as a 
non-dynamical prescribed background) or no preferred metric is available one may find 
alternative approaches leading to conserved quantities. In such circumstances there exist 
3-forms derived from a diffeomorphism invariant actions J M A m that give rise to certain 
vector valued maps on vector fields. These give rise to two other maps T N and T B which 
have traditionally been associated with "canonical stress-energy-momentum tensors" and 
are related by the so called Belinfante-Rosenfeld procedure. It is shown below that these 

3 In view of the ambitions of M-theory this interpretation may be an effective one imposed by dimen- 
sional reduction. 



3 



give rise to conserved currents in the presence of background fields which possess appro- 
priate Lie-symmetries and that T n gives rise to conserved currents, in general, only if g is 
dynamical and satisfies Einstein's theory of general relativity. Furthermore it is demon- 
strated that any difference between T B and T n arises from the dependence of A m on the 
background fields {Z}. In this article the consequences of diffeomorphism invariance of 
matter actions on all these quantities will be explored. 

The formalism below first establishes an intrinsic variational calculus for actions involv- 
ing the metric tensor field g, arbitrary tensor fields {Z} and differential forms {£}, {c?C} 
on an n-dimensional manifold M. Any differential form of degree n on M will be called 
a top-form. The set of all p-form fields is written TA P M. Thus the matter action is some 
top-form: 

A m (g, Z 1 ,...,Z N ,(i,d( 1 ,...,( Q , d( Q ) G TA n M (1) 
and depends in general on 

• a metric tensor field g, 

• a collection of mixed degree tensor fields (including scalar fields) of no particular 
symmetry Za for A = 1 , . . . , N 

• a collection of differential forms Cb G TA Pb M for B = 1, . . . , Q and their exterior 
derivatives d(s- These fields will be taken to satisfy variational field equations 
following from some action top-form A m . Since A Ein is independent of Za and £b 
these are the same as the variational equations that follow from the action top-form 
A T . 

The notation used in the paper is given in section 11.11 In section [2] partial Gateaux 
derivatives of A m with respect to tensors and differential forms are related and compared 
with standard partial variational derivatives using an intrinsic formulation. Diffeomor- 
phism invariance of A m is defined and its consequences expressed in terms of certain 
maps A and B on these derivatives. In section [3] the Einstein-Hilbert and other stress- 
energy-momentum tensors are defined and expressed in terms of these maps. Relations 
are derived between these tensors and tensor densities when all matter fields satisfy the 
variational field equations derived from A m . It is then shown how conserved quantities 
can arise in the presence of material and Killing Lie-symmetric background fields. In 
the concluding section the physical implications of these relations are emphasized. The 
Appendix gives some mathematical details and proofs of results used in the main text. 

1.1 Notation 

Local coordinates on M are denoted ( X ■ • • • • X ). These define a local coordinate frame 
{<9i, . . . , d n } where d a = and coordinate co-frame {dx 1 , . . . , dx n }. In these frames a 
metric g = g a bdx a ® dx b and the inverse metric g = g ab d a ® db where g ab gbc = &c- Here 
implicit summation is over a,b = 1, . . . , n. The metric dual of any vector field v G TTM 
is the 1-form v — g(v,—) G rA 1 ^^. In spacetime n = 4 and g is Lorentzian with signature 
(— 1, +1, +1, +1) here. A map T : X H> T(X) is /-linear if T(fX) = fT(X) for all 
scalar fields / G rA°M. The map T is IR-linear if T(XX) = XT(X) for all constants 
A G R. In the following the word "tensor" refers to an /-multilinear map on vectors and 
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their duals (co- vectors). Given a non- vanishing top- form Q, the phrase "T is a tensor 
density with respect to Q of weight W" implies that when Q is replaced by Cl = JQ for 
J E rA°M nowhere vanishing then T transforms to T = J W T. If = dx 1 A • • • A dx n 
and f2 = dy l A • • ■ A cfa/ n are related by coordinate transformation then J is the Jacobian 
of the transformation. Choosing f2 = *1 one can convert T into a bona-fide tensor field 
(T = T). The Lie derivative of any tensor T with respect to any vector field v on M is 
denoted A,T and the exterior derivative d on differential forms is defined so that d 2 = 0. 
A form (3 is said to be closed if dj3 = and exact if /3 = da for some a. The interior 
contraction operator with respect to v on forms, denoted i v , is a graded derivative and 

W 2 = 0. 

2 Intrinsic Variational Calculus 

2.1 Algebraic Preliminaries 

The degree of an arbitrary tensor will be represented as an ordered list s of or more en- 
tries. Each entry is either the symbol F (for 1-form) or V (for vector) e.g. s = [¥, V, V, F]. 
The bundle of tensors of degree s over M is denoted ® S M, with sections in T(g) s M and 
the bundles of 0-forms, 1-forms and vector fields are written 

A°M = <g) [] M , A'.U = <g) [F] M and TM = ® [V] M 

respectively. Furthermore 

(® S M) ® (®*M) = ® [s ' t] M 
where [s,t] is simply the concatenation of the two lists. Thus 

t)®(®«e r<g) [v ' F - v] M 

where v,u E TTM and ( E TA 1 M. The metric g lies in the symmetric sub-bundle of 
(g)[ F ' F ]M and the inverse metric g lies in the symmetric sub-bundle of (gP^M. Similarly 
since 

W A '" Ae ^l ^ ^a/^fjje^^-^e^ (2) 
h<...<i p P ' h<...<I p aeS p 

where 5 P is the set of permutations of {In, ... , J p } and e(cr) is the signature of the per- 
mutation, the j9-form bundle A P M is the antisymmetric sub-bundle of (gj^' - '^M where 
[F, . . . , F] has length p. 

The dual space of ® S M is ® S M where s is obtained by interchanging the symbols 
F and V in s. The total contraction of elements in T® S M with elements in T(g) s M is 
written 

T®°M x T(8) S M ->■ TA°M , ($,Z)^$:Z 
where $ G T(8)*M and Z G T® S M. It is defined inductively via 

= v: ( = (( v ) where C e T(8) [F] M and u G 1(g) .1/ 
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and extended by /-linearity to arbitrary tensors via 

($1 ® $ 2 ) : (Zi ® Z 2 ) = ($i:Z 1 )($ 2 :Z 2 ) 

Thus for example 

(Ci © u ® C2) : («i ® /3 ® u 2 ) = Ci(«i) £0) Czfo) 

where (Ci ® u © C2) e I^^M and <g> /3 ® u 2 ) G r®[ v < F ' V ]M. 

Using this notation A m may be regarded as a fibre bundle morphism 

A m : S i9) © £ {z) © £ {0 ->■ A n M (3) 

which we write A m (#, Zi, . . . , Z N , G, dCi, • • • , Cq, <Kq). Here 

• gia) = ®MAf is the bundle of metrics, i.e. G T£^\ 

• £^ = (g) Sl M © • • • © <S) SN M is the bundle of A?" tensors of the appropriate degrees, 
i.e (Z u ...,Z N )e T£^ where Z A G T(£) SA M. 

• £^ = A P1 M © A P1+1 M © • • • © A P «M © A P « +1 M is the bundle of Q pairs of forms 
of the appropriate degrees, i.e (Ci, d(i, . . . , £q, c/(q) ^ where G TA Ps M. 

Since A m is a fibre bundle morphism, the value of A m (g, Z x , . . . , Z^r, 

(i,d(i, . . . , Cq, dC,o)\ x for some point x G M depends only on the values of its arguments 

at that point, i.e. 

A m (#, Zi,..., Z N , Ci, dCi, ■ ■ ■ . Cq, dCo)U = 

^■ m (^U) ^iUj • • • ) -ZjvU, Cilacj ^CiU, ■ ■ ■ , CqIx> ^CqU) 

In this language a model containing matter in a linear (temporally and spatially) non- 
dispersive medium interacting with the electromagnetic field F = dA is characterised by 
a U(l) gauge invariant excitation tensor G = x(9j Zi, . . . , Zn, dA, (2, d( 2 , ■ ■ ■ , (q, c^Cq) f° r 
some constitutive tensor x an d is described by an action 4-form on spacetime Ai 

A m (g, Z X ,...,Z N ,A, dA, ( 2 , d( 2 , . . . , Cq, d( Q ) = 

±FA*G + A Cl (g,C2,D( 2 ,...,( Q ,D( Q ) 

where D£ = d£ for electrically neutral real fields £ and is the U(l) exterior covariant 
derivative for complex charged fields. 

A particular model [2J 13] involving only a single non-dynamic tensor Z together with 
F and g is described by the action 4-form A m : £^ © © -> A A M (with the 
bundles £^ = & ¥ > ¥ > Y WM, £^ ] = A 2 M over spacetime M): 

A m {g, Z, dA) = jF A *(Z(F) + Z\F)) (5) 

with 

F = dA (6) 
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and where is the adjoint of Z defined by 



aA^)=^A*Z(a) (7) 

for all a, (3 £ TA 2 Ai. Unlike Z, the tensor Z* depends on g. Varying (j^J) with respect to 
A then yields, in terms of the notation defined in f |22|) below 



5A m 

— (jg t Z,dA) = 

i.e. 

d*G = 

where in terms of fl20l) below, the excitation tensor G G TA 2 Ai is given by 



l.C 



G = X {F) = \(Z{F) + Z\F)) (10) 

Equations (JBJ) and (JHJ) constitute the "on-shell" Maxwell system for model in any 
background g, Z. In this model the constitutive tensor x is independent of the motion of 
the medium. Although from (J7|) Z^ depends on the metric it follows from (1T0]) that (151) 
can be written 

A m (#,Z,tL4) = |FA*Z(F) (11) 

and thus the only metric dependence of A m is through the Hodge map. 

A more complex model [3] in which the constitutive tensor x depends explicitly on 
the motion of the medium and exhibits intrinsic magneto-electric constitutive properties 
involves a timelik^] vector field V and four background degree 2 tensors Z de , Z db , Z he , Z hb . 
It is described by the action 4-form A m : E {g) © S (z) © S (c) ->■ A 4 A4 where = 
®V™M © <S>^M © ®V™M © gP^A* © g^Al, £(0 = A 2 M, 

A m (g, Z dc , Z dh , Z hc , Z hh , V, dA) 

= \F A *(z de (i Vg F) A V g + Z dh {i Vg * F) AV g (12 ) 



where F = dA, 



*{Z he (i Vg F) A V g ) - *(Z hh (i Vg *F) A V g ) 



V 9 = . y== , V g = g{V g -) (13) 

V-g(V,v) 



Since small variations in the metric g do not change the timelike nature of V, its Gateaux derivative 
with respect to g is zero. 
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and the subscripts g indicate explicit dependence on the metric. Again variation with 
respect to A gives the Maxwell equation (j£J) where G is given by flHJ). Thus 



G = Zf{i Vg F) AV g + Zf{i Vg * F) A V g 
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*{Zf{i Vg F) A V g ) - *{Zf{i Vg * F) A V g ) 



(14) 



where 



^ de = K ° (^ dc + (z dc Y) o tt, , zf = o (z db - (z hc )t) O ^ 9 , 

^ = ¥9 ° - (Z dh y) o tt 9 , Zj b = §7r g o (Z hb + (Z hb )t) o vr, , (15) 
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Id; • \; 7 \; (/ and aA^Zj(/3) =/3A*(Zj) t (a) 



for G {Z dc ,Z db ,Z bc ,Z bb } and a, (3 G TA 1 M. This implies Z\ is spatial with respect 
to V g , i.e. Z l g {V g ) = and iv g Z l g (a) = for all a G rA 1 ]^". Since vr 9 (a) AV s = aAF 5 and 
TT g (i Vg -f) = i Vg -f for all a G TA X M and 7 G TA 2 M then Z de may be replaced by \{Z dc + 
(Z dc )t) in G. Similarly Zf may be replaced by \ (Z dh - (Z hc )t), Z ff hc by \ (Z hc - (Z db )t) 
and Z bb by §(Z hb + (Z hh YY This expresses G more simply in terms of the constitutive 
tensors in the action top form (\Y2\i . Furthermore, after some rearrangement, one finds 
that \F A*G = A m , (cf [31). 

To facilitate the presentation below, it proves useful to relabel tensors in the arguments 
of A m as 

Zq = g , Z A = Z A , Z N+2B -\ = Cb and Z N+2B = d( B (16) 

for A — 1, . . . , N and B — 1, . . . , Q, so that 

A m (Z ,...Z N+2Q )eTA n M (17) 

The range C = 0, . . . , N + 2Q will be used to index the Zq- 

In the following tensors of the form * = fi<g>$ G r(A n M®® ¥ M) where fi G TA n M and 
$ G r® s M arise naturally by "differentiating" A m with respect to one of its arguments. 
One may contract such a tensor with Y G T0 s M to isolate Vt 

r(A n M ® <8) ¥ m) x r® s M -> ta"m , (^,r)^*;y 

according to the rule 

(fi®$);y = ($ : y)n (is) 

where G TA"M and $ G T® ¥ M and F G r<g) s M 
2.2 Variational Derivatives 

AA m 

Using the indexing notation (I16p and (IT7|) . the Gateaux derivative G r(A n M ® 
' ^ AZ C 

SC M) of A m with respect to Z c G r<g) Sc M is defined so that 
AA m ._ d 



■:Y 



AZ C de 



e=0 



A m (Z , . . . , Z C -i,Z c + sY, Zc+i, • • • , Z N+2Q ) (19) 
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for all Y G T(g) Sc M. An example in a local frame is given in appendix [B] 

By contrast, for C = N + 2B - 1 then Z N+2B -i = (b G TA Pb M i.e. Cb e r®^--^ 

<9A m 

(with a list of length p B ) the Gateaux derivative — — G TA n Ps M is defined so that 



a A 



<9A r 



<9Cb de 



e=0 



A m (g, Zi,..., Z Nl Ci, d(t, . . . , ( B + ea, d( B , ■ ■ ■ , (q, d( Q ) 



AA r 



(20) 



AC 



-:oc 



B 



dA 



for all a G rA Ps M. The derivative — — G TA n PB M may be related to the derivative 



AA r 
AO 



dc 



B 



g r(A n M 



since one can identify TA n ~ PB M and {fi g> V G T(A n M 



|V antisymmetric}. The correspondence follows from the relation: 

(fi® V):a = a Ai v fi (21) 
for any a G TA P M. Here a general antisymmetric tensor V G r®^ v ' " ,V ^M can be written 

V = ^ V Wp 4A.-Al; p 

h<...<I p 



(22) 



pi 



Il<...<Ip CTgSp 

and the internal contraction operator iy with respect to V is defined by: 



ii<...<J p 



so that for a G TA P M 



a: V = iv tt 



(23) 



(24) 



The proofs of (j2ip and (1241) are given in lemmas [7] and [6] respectively in appendix lAl 

<9A m 

Likewise if C = N + 2B then Zjv+2B = d( B G rA Ps+1 M and tttt^ G rA n_Pfl_1 M is 



defined by 



d(d( 



B, 



/3 A 



dA r 



a(dc 



B 



de 



£=0 



A m (g, Z 1 ,...,Z N ,( u d( 1 ,...,( B , d( B + e(3,..., ( Q , d( Q ) 



AA 1 



(25) 



A(d( B , 
for all G rA Pfl+1 M. 



■0 



5A 1 



The variational derivative — — of A m with respect to Cb G rA Ps M is defined by 

<Kb 



„ 5A m d 
M <Kb de 



6 = 



A m ( 5 ,Z 1 ,...,Z iV ,Ci,dCi,..., 

Cb + c^Cb + • • • , Cq, ^Cq) 



(26) 



for all j3 G rA Ps M with compact support. Hence 

SA m dA m 



5Cb d( B \d(d( B 

follows from lemma [8] in appendix |A] 

A p-form (b is said to be "on A m -shell" if 

Kb 

In this situation 



+ ( - 1)PS+1 <a^r) (27) 



:28i 



<9A m / <9A m \ 

dCs ^d(dCij)/ 
The matter system is said to be on A m -shell if (T2"gj) is true for all 

2.3 Diffeomorphism invar iance 

Given a local diffeomorphism (ft : t/j^ — > U' M where Um, U' m C M then the two maps 
4>± and induce the map : <S> s ^m — > ® S U' M on tensors. The action top-form 

A m (Zo, . . . Zn+2q) is (locally) diffeomorphism invariant if 

^(A m (Z , Za, . . . , Z N+ 2 Q )) = A m (0(Z o ), i(Z x ), ^(Z N+2Q )) (30) 

for all local diffeomorphism <fi : Um —> U' M . 

Lemma 1. If A m is diffeomorphism invariant then for all vector fields v G TTM 

N+2Q AA m 

£„(A m (Z , Zi, . . . , Z n+2 q)) = ^2 ~\Z~'^ vZc ( 31 ) 

c=o 

which may be written 

C v (A m (g, Z 1 ,...,Z N ,(i,d(i,---,(Q, d(g)) = 

AA m . AAA m .„ ^ „ <9A m ^, Mf , 5A m (32) 

_ :£ „ s + g _ :£ . Z , + g « B A _ + g £ „ (dCB ) A ^ 

Proof. Let e be the one parameter family of diffeomorphisms generated by v. From (1301) 

d 



C v A m (Z , . . . , Z n+2 q) — — J) e (^A m (Z , . . . , Z 



d_ 

de lo 

N+2Q 



N+2Q, 
e{^0), ■ ■ ■ , 4>e(2 N+2 Q i 



— A m i<bJZ n ),...,6JZ 



E "^| o A m (^,...,0 £ (2c),...,2 



c=o 



JV+2Q 



E X (^°' • ■ ■ ) -^c + eL v Zc, ...,Z_ 



c=o 

N+2Q 



N+2Q 



x ^ AA m 



c=o 

Then (El follows from (US) and (ESI). □ 
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<9A m <9A m 
Since £ v (b = MCs + di v ( B , the top-forms £ V ( B A ^— and C v {dC, B ) A . - . in (SB 

can be expressed as 

<9A m / <9A m \ <9A m <9A m 

£„Cb A — - = d (i v ( B A — — + (-1) Ps ^Cb A d-— + i v d( B A — - (33) 
dC B \ d( B ) d( B d( B 

and 

dA m / dA m \ dA m 

c " {d(B) A WdcT) = T A A awc^) " ( - 1)M '" dCB A (34) 



where Cb G rA Ps M. See lemma [9] in the appendix lAl 

AA m . 

By analogy with ( !3"3| and (]3"1|). the top-forms -C v Zc in ( T3~Tj) may be written 



AA m . 



/ /AA m \\ /AA m 



(35) 



AZ C 
where the maps 

A : TTM x T(A n M ® ®*M) x T(g) s M -> rA n_1 M , (v, Z) i-> ^(^, Z) 

and 

B : YTM x T(A n M <g> ®°M) x r® s Af -»■ TA n M , (v, Z) i— >■ , Z) 
are 'f'-linear in i> and defined inductively as follows: For 0-forms, / G r<g)^M 

A(^,/) = and /) = (36) 

For 1-forms aeT®% 



® w, a) = a(v)i u Q 
and H„(fi ® u, a) — (—l) n+1 i u Q A i v da — di u Q A z„Q! 



(37) 



For vectors u G T(g) [vl M 



A v (fl <S> a,u) — —a(v)i u Q 
and B v (Vt ® a, u) = v(a(it)) f2 — A i^c/a + (— l) n di u f2 A 2„a 



(3* 



For tensors $i G T(g) s M, $ 2 G r®*M, Zx G r<g) s M, Z 2 G T0*M and tt G rA"M, define 
the Leibnitz rules: 

Av{n ® $i ® $ 2 , Zi <g> Z 2 ) 

= A(($i : ® $2, Z 2 ) + A(($2 : Z 2 )ft ® $i, ^i) 

fi„(f2(8) $i (g) $ 2 ,Z! <g> Z 2 ) 

= : <g> $ 2 , Z 2 ) + ^(($ 2 : z 2 )n ® $x, zo 

The proof that (|35j) follows from lemma [TD] given in appendix E The proof that (133 p and 
AM]) are consistent with ( 135]) is given in lemma [IT] appendix |Aj 
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If A m is diffeomorphism invariant then from fl3T|) and f |35|) : 

N+2Q AA m 
^A m = £,A m = V -— - : C V Z C 

c — o c 

= 4 +E*b* 



( '-(> ° 7 c*=o 

hence applying lemma [5J appendix [A] gives 



and 



and 



*(^.')-f:^^)-i:((-^^ 

y A=l B=l v ^ 



0Cb 

m 



<9A m dA 

+ i v d( B A -5— + (-l) Pfl+1 ^dCs Ad 



5Cb v y b " 9(dC B : 

When the system is on A m -shell (I43p may be written 



, /AA m \ tm , , 5A m 



. . dA m \ A . /AA m _ 

whereas (jUJ) reduces to 



a=i — ' A 



(40) 



z,A m = £ A(^Uc) (41) 



Using (1551) and (}3"4"|) . these may be written as 

/AA m \ N /AA m 

' A=l A 

<9A m _ <9A 



(43) 



(44) 



(45) 



See lemma [T21 appendix IA1 Relations (|43|) - (|46p play a pivotal role in the arguments below. 
In particular they enable one identify terms which may be associated with the quantities 
derived historically with Noether, Belinfante and Rosenfeld. 
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3 Currents and Conservation laws 



Many quantities in physics owe their raison-d'etre to the existence of conserved quantities 
that do not change with time in a dynamical system. Thus notions of energy, momentum 
and angular momentum arose from the analysis of Newtonian particle dynamics. With 
the introduction of fields and the development of continuum mechanics it became natural 
to incorporate such concepts into continuous dynamical systems and their unification into 
a "stress-energy-momentum" complex offered an attractive objective. However, as is well 
known such a unification is not unique and the indiscriminate use of the term "stress- 
energy-momentum" tensor has led to unnecessary confusion when discussing forces and 
torques produced by fields in media. 

In this section a number of technical issues are addressed that inter-relate these physi- 
cal concepts. They include the role played by different aspects of (multi-)linearity needed 
for a general definition of stress over curved surfaces in space, the role of Killing symmetry 
needed to establish the notions of energy and momentum and the role of algebraic sym- 
metry of maps and associated tensors or tensor-densities in their conservation. Using the 
variational framework established in the previous section it is possible to correlate these 
aspects with the parts played by the presence or absence of background matter fields and 
background gravitation in their implementation. 

Since some of the traditional arguments for the construction of a symmetric stress- 
energy-momentum tensor in Minkowski spacetime are spurious (even in the absence of 
background matter fields) it is useful to begin the discussion with the Noether current 
associated with A m in a general background using the Lie- derivative. This leads naturally 
to conservation laws in the presence of background symmetries. In this manner it is also 
straightforward to extricate the role played by /-linearity in establishing a proper tensor 
description of stress. 



3.1 Noether and Belinfante-Rosenfeld currents 

For the restricted class of actions considered in this article we define the Noether (n — 1)- 
form current associated with A m by 

K = i v A m - V £ v (b A — — (47) 

for v G TTM, where it is assumed that the system is on A m -shell. Since the Lie derivative 
C v is not /-linear in v neither is M v and for v = v a d a , lemma [TBI in appendix lAl gives 

dA m 

K = v a M 9a + Y,dv a A iauCfl A j^— (48) 

It is however M-linear in v (where v a are constants). For a chosen nowhere vanishing 
Q E TA n M one may define the map T N : TTM -»■ TTM by 

Af v = i T M {v) VL (49) 

In a coordinate system (x 1 , . . . , x n ) with f2 = dx 1 A ■ ■ ■ A dx n and A m = SVL then has 
component maps 

T Ma h = dx a (T M {d b )) (50) 
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which from lemma [15] in appendix [A] gives 

dx a AAf db = T Na b n and T Ma b = i dn ---i dl (dx a AAf db ) (51) 
Lemma [27] in appendix [A] then gives 

- <9<? 



where 



b=i h<-<i PB 



*, B = ^2 Cbi dx h A • • • A dx Ip s 
h<-<i PB 



bmce T N is not /-linear, i.e. T N {fv) ^ fT M ( v) in general, one must not confuse b a 
with the components of a tensor field. In fact using (T4"gj) and lemma [TS] in the appendix 
for v = v a d a , 



Q 



T N (v) = (T Af ) a b v b + (d c v b )i dl ■■■i 9 Jj2 dx a A dx c A i db C B A — — ) (53) 

0{dC, B )/ 

However, it is not uncommon to refer to ( 15"2"|) as the components of the canonical stress- 
energy-momentum tensor associated with A m in Minkowski spacetime [U [HJ [13] . This 

x a -> y a = R a b x b + 5 Q 



arises since l~^b a R a cT^b c , under affine coordinate transformations of the form 



where det(i?) = 1 and R a b, are constants. This accounts for its widespread use in 

special relativity as formulated by Einstein and Minkowski in spacetime Ai. Since is 
not tensorial with respect to arbitrary coordinate transformations its use for calculating 
stresses is restricted to planar surfaces in space. This follows from (1531) . since for any event 
p E Ai on a non-planar spacelike 2-surface with normal field w, T^(w)\ p will depend on 
the derivatives (d a w b )\ p . Thus one cannot, in general, define Cauchy traction forces that 
must be independent of such derivatives [12] . 

The requirement that the concept of stress follows from a bona-fide tensor leads one 
to seek an /-linear map constructed from M v . Since C v = di v + i v d, the (n — l)-form 
current 



r. 



f = K + dS v (54) 



where 



B=l 



is manifestly /-linear in v. One may refer to this as a Belinfante-Rosenfeld formula [5] 
although one may note that no metric on M is necessary for its construction. It follows 
from gTJ) and <^ that 

r° = a- - E (w. a + wc» A 5(35)) (56) 
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which we refer to as the Belinf ante- Rosenf eld stress-energy-momentum current associated 
with A m . From f )43|) one can also write fl56j) in terms of the maps A v as 



A=l 



(57) 



Again, for a chosen nowhere vanishing VL e TA n M the map 

r e e r ®[v,F] defineci by 



«r s (u) 



(5f 



is a density with respect to with weight —1. 

Given a preferred metric g and any map T : TTM — > TTM one may define the map 
T : TTM x TTM TA°M by 

T(u, v) = g(T(u),v) 

This enables one to discuss the algebraic symmetries of T. Such a map is said to be 
algebraically symmetric with respect to g if 



T(u, v) = T(v, u) 
This implies (see lemma [TH in appendix [AT) 

v A ij(u)Q — u A ir(v)Q 



(59) 



(60) 



for all it, f G TTM any non- vanishing top-form G A n M. Thus for the maps T - ^ 
FTM x TTM -»■ TA°M and T B : TTM x TTM TA M 



T^(u,<u) = ^(T^(u),i;) and T B (u,v) = g(T B (u),v) 



Then T-^ is symmetric if 



u A JV« = v A JV« for all FTM 



and T is symmetric if 



u A -r 



B 



v A r,. 



for all 



u, v e TTM 



(61) 



(62) 



(63) 



To illustrate these notions, consider the premetric formulation of electromagnetism 
0, [6] on spacetime M. where one starts with the action 

A m (Zi, . . . , Z N , dA) = \F /\H (64) 

with F = dA and % = k(Zi, . . . , Zn, F) linear in F. The Noether current is then 

N v — \Ff\ i v H - \i v F AU + di v A A U (65) 



5 Some authors [111 113) refer to T as the canonical stress-energy- momentum tensor whereas others 
[TUJ [7] refer to T B as the canonical stress-energy-momentum tensor. The appellations T 1 ^ and T 8 
eliminate this notational ambiguity. 
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This is manifestly not U(l) gauge invariant. However the Belinfante-Rosenfeld current 



T. 



B _ 1 



F A i v % — \i v F A U (66) 



is U{1) gauge invariant. 

When a metric is prescribed as in the spacetime model (jl]) where 

A m (g, Z u ..., Z N , dA) — \F/\-kG (67) 

with G = x(9i Zi, . . . , Z N , F) linear in F, the Noether current is given by ( 16"51) with 
H = *G and the Belinfante-Rosenfeld current is given by ( 166|) with H = *G. In this case 
the Belinfante-Rosenfeld current 

T v = If A i v * G - iz„F A *G (68) 

gives rise via fl58l) and fl59|) to the algebraically non-symmetric Minkowski stress-energy- 
momentum tensor [9] density T B . Examples (jSJ) and f|T2|) are particular cases of fl67|) . 

In general, neither T-^ nor T B possess the algebraic symmetry (|59|) due to the presence 
of background fields. In the particular case of the vacuum where % = *F, then T B has 
algebraic symmetry, while remains non-symmetric. In the absence of a preferred 
metric one cannot even define or T B from and T B respectively. 

3.2 Conservation Laws 

In the presence of Lie-symmetries of g and Za, both the Noether and Belinfante-Rosenfeld 
stress-energy-momentum current give rise to conserved material quantities. 

Theorem 2. If K G YTM is a Killing vector field, i.e. Cxg = 0, and in addition 
the background tensor fields {Z} satisfy the Lie-symmetry condition CrZa = for A = 
1, . . . , N, then both the Noether current and the Belinfante-Rosenfeld current are closed 

d T % = dU K = (69) 

Proof. Since C K g = and CrZa = then from (J35J) 

/ /AA m \\ /AA m \ 
d iAK[^r,g) ) +B K [^r,g)=0 (70) 

and 



V Ag >») ' ~ A \ Ag 



( /AA m \\ /AA m \ 

%Ma^> z -))+Maz?^) 



for all Za- Thus from (1571) and (146 



(/AA m \\ N ( /AA m 
H^ 9 )) + ^{HaZ-a> Za 

f AA m n A / AA m 



A=l 



□ 
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In terms of the maps and T B , using lemma [TBI in appendix |A]> (1691) may be written 
in a general coordinate basis with Q = dx 1 A ■ ■ • A dx n as 



d a {{T B ) a b K h ) = (71) 



and from (153 



/ Q 
d a I {T M )\K b + {d c K b )i dl ... ldn ^dx a A dx c A i 9b ( B A 

^ B=l 



aA ' - (72, 



d(dc 



Bj 



When g is Lorentzian on spacetime, by Stoke's theorem both T B and yield conser- 
vation laws. An energy conservation law follows if K is timelike, a linear momentum 
conservation law if there exist three independent spacelike translational Killing vectors 
and an angular momentum conservation law follows if there exists a basis of three in- 
dependent spacelike Killing vectors generating spatial rotations. Such conservation laws 
make no reference to forces (stress) or torques (moments) and are valid in the presence of 
smooth non-dynamical background fields. However only T e , being /-linear, deserves the 
appellation stress-energy-momentum tensor. 

If A m does not depend on a metric g the requirement that K is a Killing vector may 
be dropped, i.e. if A m (Z 1 , . . . , Z N , d, d(i, . . . , (q, d(g) is independent of g and CyZ^ = 
for A = 1, N and V e TTM, then 

dr® = dM v = (73) 

However, in the absence of a metric, no physical concept of energy or momentum exists. 

If Ck9 = but not all background fields {Z} are Lie-symmetric then a simple gener- 
alisation of theorem [2] yields 

(/AA m \ /AA m 
dAx ( A2~' Za ) + Bk \~KZ~' Za 



{A\C K Z A ^0} 

AA m . 

AZ A 



E w?: :XkZa 



(74) 



{A\C K Z A ^0} 

which in general is not equal to zero. 



3.3 Historical perspectives in Minkowski spacetime 

In the absence of background fields {Z} in Minkowski spacetime, T B ^ = T B ba and from 
f lTT]) . d a {T Ba b) = 0. It follows that the "moment of T v q" is conserved. Thus in inertial 
coordinates ( 3 ) with 

M n U = —(x^T^vq — x v T v q) 

for fj,, v — 1, 2, 3, vanishing at spatial infinity, one has from f l7Tj) with K = x^dy — x u d^ 



d 
dx° 



/ M E '^dx 1 A dx 2 A dx 3 = (75) 
Jr 3 
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This is identified with the conservation of (orbital) angular momentum of a field system 
in M 3 . This result is a direct consequence of flTTj) which does not require the imposition 
of any algebraical symmetry 

More generally with the Killing vector fields, for a, b = 0, 1, 2, 3 

Tlab = x a d b - x b d a (76) 

and 

T B {R ab ) = M B ab c d c (77) 

where 

M B ab c = T Bc a x b - T Bc b x a (78) 
is the "moment of T Bc b , one has 

d c M B ab c = (79) 

In fact (ITTj) shows that f[77)j) and f[T9"j) are valid in the presence of Lie-symmetric background 
fields satisfying L-ji^Za = for all A and for all 1Z ab . 

In Minkowski spacetime there exists a basis of 10 Killing vector fields generating the 
algebra of the Poincare group. The nature of the algebraic symmetry of T B ab depends on 
the Lie symmetry of the background fields. This is a consequence of the following lemma: 

Lemma 3. In Minkowski spacetime, for fixed a and b, let d a and d b be two commuting 
translational Killing vectors such that Cg a ZA = and CQ b Z^ = for all A. Then T B is 
algebraically partially symmetric if and only if T~n ab is closed, i.e. 

T B ab = T B ba & ^L = ° ( 8 °) 

Proof. Since CQ a Z^ = and Cq^a = then drg a = and drg b = 0. From /-linearity 
T R a t = XaT d b ~ x b T d a h ence from (I6"0"|) 

dr% ab = dx a A r| - dx b A r| o = {T B b - T b B )dx° A dx 1 A dx 2 A dx 3 
Hence (EDI. □ 



Thus total (orbital) angular momentum conservation in Minkowski spacetime, gener- 
ated by {di,d 2 ,d 3 }, does not demand that T B ab is fully algebraically symmetric for all 
a, b. 

It is of interest to note that conservation of an 5*0(3, 1) Killing current can be obtained 
directly from the Noether quantities in the absence of background fields. It follows from 
(l69~j) and flTT]) that for all Killing vector fields lZ ab one has 



dMn ah = (81) 

and 

d c {T M {R ab ) c ) = (82) 
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where from (j48l) 



^ / dA m dA 

Nn ab = x a H db - x b Af 9b + ^ [ dx a A id b ( B A - dx 6 A id a (,B A 



m 



0(dC B ) ^ 5(dCi 

and from 

T N (R a b) = T jyc a Xb — T Nc h x a 

+ cfe c A ete a A z^Cb A — - - dx c A dx b A i 5a Cs A — — 

Thus the Noether current A/7? , does not coincide with the moment of Noether linear 
momentum: x a Afg b — Xf ) Afg b . The additional terms result from the fact that M v is not f- 
linear in v. However although both T B and give rise to conserved quantities, the lack 
of /-linearity in precludes its use for the definition of stress over curved 2-surfaces. 

Equations (1HT1) and (|82|) remain valid in the presence of Lie-symmetric background 
fields with C-ji^Za = for A = 1, . . . , N and fixed a, b, 

3.4 The Einstein-Hilbert stress-energy-momentum tensor and 
its associated currents 

AA m 

In general relativity the variational derivative — — is used to define the algebraically 

Ag 

symmetric Einstein-Hilbert stress-energy-momentum tensor T H E T^ F,¥ ^M for the the- 
ory: 

/ AA™ \ 

T H {u,v) = 2*- 1 Ku®v)) for u,veTTM (83) 

which is manifestly /-linear in u and v. 

Using an arbitrary vector field v on M with a metric g, it is convenient to use T H to 
define the associated (n — l)-form current G rA n_1 M by 

7* = *(?>,-)) (84) 

which is manifestly /-linear in v. Thus 

T n = (*-V«) 9 e a (85) 

where {e a } and {X a } constitute mutually dual frames. From ( 183]) it is also clear that 
since g is a symmetric tensor, T H satisfy the algebraic symmetry 

T n (u, v) = T n (v, u) for all u,v E TTM (86) 

(c.f. fl59|) ) and hence {rjf } satisfy the algebraic symmetry condition 

e a A - e b A r£ = (87) 

where e a = g(X a ,X b )e b . 
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In lemma [16] appendix [A] it is shown that 



/ AA m 

* = Aj— — ,g) for all v e TTM 



Thus from (143]) it follows that in terms of the map A v , 

N /AA m \ / <9A m <9A m \ 

r« = ,„A™ Zx) - £ («. A ^- + Wfe A — ) (89) 

Let the exterior covariant derivative of be defined by 

(Dt h ) v = {i v e a )dr^ a - i v de a A r« (90) 

This corresponds to the standard Levi-Civita covariant derivative in the case of a torsion 
free metric-compatible connection V and (Dr n ) v = implies V • T n = 0. 
In lemma [T7I appendix ]A1 it is shown that in terms of the map B v , 

(DT n ) v = -B v (—,g) (91) 

Hence for a system on A m -shell one has, from 

N 



(^ w ),=£^(=£r-,^) (92) 



A=l 



It is worth stressing that, in general, even if N = this relation is not a conservation 
law since in general (Dt h ) v ^ dr£. However, if K e TTM is a Killing vector field, 
C K g = 0, then from (JHHJ), (]&]) and (]9B 



i.e. 

(£>r*)* = dr£ (93) 

From the definitions above, the relationship between the Einstein-Hilbert, Noether 
and Belinfante-Rosenfeld currents is given by 

N /AA m \ 
r« + X) A f^- ; ^a) = = jV« + dS, (94) 

A=l A 

As stated in (1691) rf and jVr- give rise to conserved currents associated with each 
Killing vector field K in the presence of the Lie- symmetric background fields {Z}. By 
contrast, from (193]) and (I9T1) . in the presence of arbitrary background fields, whether 
Lie-symmetric or not, 

dr% = f:B K (^f,Z A ) (95) 

A=l A 



20 



which in general does not vanish. However in general, as stated above, unlike neither 
Tjf nor Mk possesses the algebraic symmetry condition f l59|) .f lH7|) in the presence of any 
background fields. 

If there are no background tensor fields {Z}, i.e. N = then from (19"1"|) 

t* = T ?=N v + dS v for all veYTM (96) 

and hence if K is Killing from f[6"9"j) 

d T % = dr® = dM K = (97) 



In this case, since satisfies the algebraically symmetry condition it follows from 
that also does. However in general M v does not. 

In the context of model (fT2l) the Einstein-Hilbert stress-energy-momentum tensor T n 

_ AA m 
is given by (1831) by evaluating the derivative — — at 

[g, Z dc , Z db , Z hc , Z hh , V, dA) = (g , Y dc , Y dh , Y hc , Y hh , W, dA) (98) 

where go is an arbitrary background metric with associated Hodge map *o an d the 4- 
velocity of the medium W satisfies g (W, W) = -1. The Y l e {y de , F db , F he , Y hb } are 
spatial with respect to W, i.e. Y 1 (^g (W, — )) = and iwY l (a) = and satisfy the adjoint 
properties with respect to g : 

a A * Y dc ([3) = A * ^ de (a) , a A * Y hh (l3) = A ★ * rKb (a) 

and a A * ^ db (^) = — /3 A * ^ hc («) ' '''' ' 

for all a, /3 G rA x M. The tensor T w is the Abraham stress-energy-momentum tensor for 
electromagnetic fields in (moving) media [21 [3] . 



4 Caveats and Conclusions 

The results in this article have direct relevance to the construction of variational formula- 
tions of field systems in spacetime where complex microscopic interactions are represented 
by phenomenological macroscopic constitutive relations between the dynamic fields {£}. 
It is tacitly assumed that the variational field equations based on the matter action J M A m 
admit such {£} as non-trivial physically acceptable solutions. This condition may impose 
constraints on the constitutive modelling. For example if 

A m (#, J, A) = \dAf\*dA + AAJ 

where A e TA 1 ^ is dynamical and J G TA 3 A4 is a background field then the A m matter 
SA m 

shell condition ^ = implies dJ = 0. The physical interpretation of such models for 

continuous media can then be facilitated by finding conserved currents associated with 
continuous symmetries. 

For the class of theories described by a diffeomorphism invariant matter action J M A m , 
where A m in general depends on a background gravitation field, a collection of background 
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tensors {Z} and matter fields {£} the differential forms r^ 1 , t b and M v are related by 
f[94"j) when all matter fields {£} are "on A m -shell" and from fl95l) . (ir^ 7^ in general. 
However, in Einstein's theory 

at = n * 1 + A m ( Sj z 1} . . . , z N , Ci, dCi, • • • , Cq, ^Cq) (100) 

where 7?. is the curvature scalar derived from g. If the metric g is dynamical as well as 
the matter fields {£} one has the additional on-shell Einstein equation 

Qv = (101) 

where the Einstein (n — l)-form Q v G YK n ~ l M for v G TTM is pQ 

g v = \R{X C , X d , X b , e a )e c A e d A z v * (e a A e 6 ) (102) 

Since the Bianchi identity for Q v yields {DQ) V = for all v G TTM, ( 11011) implies 

(Dr n ) v = (103) 

Then from (l§3"j) . if AT is Killing one has the conservation law 

dr% = (104) 

and hence from (l§2"j) and f 1 1 3 j) 

A„ /AA m \ 

£Maz? z *)=° 

A=l 

This relation may impose constraints on the dependence of A m on the non-gravitational 
background structure {Z}. For example in the analysis of a gravitational wave propagat- 
ing in a material medium. 

In physical applications it is often convenient to break a system into weakly interacting 
subsystems in order to analyse the dynamics of subsystems perturbatively. However care 
is required in extending the consequences of diffeomorphism invariance for A m to top- 
forms A s describing diffeomorphic invariant sub-systems. For each A s one may define an 
associated Noether form J\f£: 

Af: = t v A s - £v<b A (105) 

Thus 

s 

Similarly 

' H -zZ<' n (107) 
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where 



/ A A s 

T s ' n (u,v) = 2*" 1 [—':(u®v) ) for a. r e VTM (108) 



and 



When all matter fields {£} are "on A m -shell" they will not necessarily be "on A s -shell"and 
one cannot associate conserved currents with subsystems, in general. This is simply 
a reflection of the interaction between subsystems in situations where g and {Z} are 
background fields. However, given some decomposition A m = X] s ^- S 5 suppose that for 

5A S ° 

some s = So there exists a subset of matter fields {Cb} such that — — = for B = 

Kb 

1 . . . Q . In this situation one may regard the quantities Cb, ^Cs ior B = Q + I . . . Q &s 

background fields to supplement those in Z. The analysis in this paper is then applicable 

by disregarding all A s with s ^ sq , replacing A m by A s ° and disregarding the equations 
5A S0 

__ ^ for B = Q + 1 . . . Q. 

Kb 

Maps such as T-^, T B and T n have been traditionally used to construct densities of 
field energy and linear momentum in Minkowski spacetime. Applied to closed systems in 
the presence of the Lie-symmetric background fields {g, Z} this procedure is strictly only 
meaningful for and T B since they alone give rise to conserved currents in Minkowski 
spacetime. Furthermore only T can meaningfully be used to describe stress and angular 
momentum in Minkowski spacetime. In a general spacetime with a non-dynamic metric 
neither T B , T n nor T give rise to conserved Killing currents unless Cxg = and 
C K Z A = for all A. 

In any background metric, the model described by (jSJ) yields a particular T B known as 
the non-symmetric Minkowski stress-energy-momentum tensor and T n as its symmetrised 
version. 

The model (Tl2|) yields a symmetric T H known as the Abraham stress-energy-momentum 
tensor (which for a general Lie-symmetric Z does not generate conserved Killing currents) . 
This model does however yield T B that coincides with the non-symmetric Minkowski 
stress-energy-momentum tensor. This sheds light on the relationship between the Abra- 
ham and non-symmetric Minkowski stress-energy-momentum tensors. The precise rela- 
tion between these two tensors, in this model follows from (1941) as 



(110) 



+ A ( +A( AAm Z^\ +A ( AAm Z^\ 

where all tensors are evaluated at the point 

{g, Z dc , Z db , Z hc , Z hb , V, dA) = (g , Y dc , Y db , Y hc , Y hh , W, dA) 

Conclusions drawn from different models are directly related to the epistemology used 
to describe the linear and angular momentum of light in unbounded media described by 
background constitutive tensor fields in background gravitational fields. 
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Thus the physical consequences of any model based on a diffeomorphism invariant 
action used to describe the dynamics of matter fields in the presence of a specified self- 
consistent non-dynamic background depend not only on the action for the model but 
also on a choice of objects, such as T H , T B or , that are adopted to define conserved 
quantities, including field energy, momentum and angular momentum. Furthermore these 
may only acquire physical cogency in the presence of sufficient Lie symmetry. 
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A Mathematical Details of results used in the text 

Lemma 4. On an n dimensional manifold M, given the (co-vector valued) forms a w G 
rA™ _1 M and j3 w G TA n M which are 'f '-linear in w such that 

d(a w ) = (3 w (111) 

for all w G TTM with compact support then a w = and f3 w = 0. 

Proof. Let w = w a X a , a a = a>x a and f5 a = /3x a - Given any subset U C M with boundary 
dU 



L*(w a a a ) = / w a (3 a 

dU JU 

where i : dU — > M is the embedding. Assume first that w has support away from the 
boundary dU and in an arbitrary small region then 

w a (3 a = 

u 

and hence j3 a = 0. Thus 

/ U/V(Oa) = 
JdU 

For all subsets U C M. Considering w to have a support on a small set about a point 
in the boundary implies t*a a = 0. Since we can choose any dU we show that all the 
components of a a = and hence a a = 0. □ 

Lemma 5. (HOD implies (jUJ and (J4"2l) . 

Proof. Since ()40p is true for all v G TTM it is true for all v = w G TTM with compact 
support. Thus ( I4ip and ( I42p follow from lemma S] setting 



N+2Q / AA m x ™ /AA 



c=o c c=o c 

□ 

Lemma 6. For any a G TA P M we have 
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Proof. Setting 



while 



a = ^2 a h ... Ip e h A • ■ ■ A e Ip 
h<...<i p 



From (ED and (E2D we have 



I\<...<Ip <JdS p 

■(h m E v^^6( P )x p(Jl) ®.--®x KJp) 

Ji<...<J p p&Sp 

w E E EE^^c-^ 

V ' h<...<I p Ji<...<J p aGSp p€S p 

h E E «r,-+V*-*4Ji-iS 
^ h<...<lp J 1 <...<jp 



h E a w„ 



p! 

^ /i<...</ p 



^« = 3 E ^ ■ ^•••^ Jl ( E a h-i p e h A • • • A e 1 * 

V' J\<...<Jp h<...<I p 

4 E E <xii...i p V Jl --- Jp ix Jp ---ix Ji e Il ^---Ae I > 

Ji<...<J p h<...<I p 



J E E aw.V*-*«g-^ 



Jx<...<J p h<...<I p 

h E <*-*v-* 



Lemma 7. For any a G TA P M, one has (l2Tj) 

Proof. For < g < p we have since • ■ • i Vq a) A (Wfi • " • V,^) e rA n M 

= (~l) q ~ 1 i Vq (i Vl ---i Vq _ 1 a) A (i Vq+1 ■ ■ ■ i Vp Q) 

= (-iy-\-i)(v-«)(i vl . . . ivq _ ia ) A i Vg (i Vq+1 ■ ■ ■ i Vp n) 

= • • ■ i Vq ^a) A (i Vq ■ --i^Q) 
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Hence from (JTHJ) and ([21 

(ft® V):a = (V:a)fi = (z v a)fi = - ^ V Wp (z Xjjl • • • i^a) A 

p " h< -<i P 

t<-<Jp 



^ h< - <l P 

a A 

F h< - <lp 

= a A Zyfi 

hence (J2TJ) . □ 
Lemma 8. (12?]) /ioWs. 
Proo/. From ([20]) and ([23]) 
oA m 



P A 

/ A m (a, Z l5 . . . , Zjv, Ci, dCi, ■ ■ • , Cb + e{3, d( B + ed(3, ...,( Q , d( Q ) 

J M 



d 



de 



, dA m IO <9A r 

/? A — — + dp A 



„ <9A m / dA m \ . „ ,/<9A m \ 

aA m 



P A 



since /? has compact support. Since this is true for all P then (127)1 follows. □ 
Lemma 9. For any £u>o forms a G rA n_p M and /3 G TA P M then 

C v p A a = d(i v p A a) + (-l) p i v p A da + i v dp A a (112) 

Proof. 

C v p A a = dzt,/3 A a + z^d/3 A a = d{i v p A a) + (— l) p i v P A da + i v dp A a 

□ 

Lemma 10. For any tensors Z G T(g) s M and * G r(A n M ® (g) ¥ M) taen 

= d(A v (V,Z)) +B V (V,Z) (113) 
Furthermore this decomposition is unique in the sense that if 

V:(C V Z) = da v + p v (114) 
where a v and P v are 'f -linear in v then 

a v = A v (^,Z) and v =B v (V,Z) (115) 
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Proof. First using ([36]) we prove that (ITT31 for / G T(g) M 

n : (A/) = n v(f) = d(A v (n, /)) + B v (n, /) 

Using (1371) this is true for 1-forms a G r®' F 'M 

(fl <g> w) : (£„a) = (i u £vu) = (-l) n+1 i u £l A c^a + A ^c?a 

= d{i u VL A z„a) — di u VL A Zt,a + (— l) n+1 i u Vt A z„<ia! 

= ® u, «)) + ® m, a) 

Now using (I38p we show that for vectors u G T^^M 

(O ® a) : (C v u) = £1 (a(C v u)) = Vl (y(a(u))) — [i u C v a) 

= f2 (y(a(u))) — d(i u VL A i v a) + di u VL A i v a — (— l) n+1 i u $l A i v da 
= d(A v (Q <g> a, «)) +B 8 (0® a, tt) 

Using fl39|) it follows that (1113P is true across tensor products. Assuming it true for Z\, Z 2 
then 

{n®^ 1 ®^ 2 ):(C v {Z 1 ®Z 2 )) 

= (n <g> $! <g> $ 2 ) : ® Z 2 ) + (fi ® $1 <g> $ 2 ) : (Zi ® £„Z 2 ) 

= : A,Zi)($ 2 : Z 2 ) + : Z x )($ 2 : £„Z 2 ) 

= ($ 2 : Z 2 )(Sl <g> $1) : (A,Z X ) + ($1 : ^i)(^ ® $2) : (A,Z 2 ) 
= d(A v {{$ 2 :Z 2 )n®§ u Z l ))+B v {{® 2 .Z 2 )n®® 1 ,Z 1 ) 

+ : ® $ 2 , Z 2 )) + : Z X )Q ® $ 2 , Z 2 ) 

= (2(A(^ ® $1 <8> $2, Zi ® Z 2 )) + B v {n ® $! (8) $ 2 , Zi ® Z 2 ) 

Hence (11131) is true for all tensors Z. 

If (11141) is true then from d 1 1 3 [) one has 

d(M*, z ) -«,) = - Z) - p v ) 

and hence from lemma H] f 1 1 1 5 j) follows. □ 

Lemma 11. Given f2(g>V G T (A n M ®<g}^ Y M) where V is antisymmetric and [V, . . . , V] 
has length p and f3 G TA P M then 

Av(P, fi ® V) = i v p A i v Q (116) 

and 

B v (/3, Vl®V) = {-lfi v p A di v Q, + i v d/3 A iyft (117) 
Proof. From ffTTBl . fl2T|) and ffTT2]) we have 

n ® V)) + 5,(0, Q <g> V) = (fi ® V) :(£,£) = (A,/?) A ? v ft 

= d^/S A 2 V ^) + (-l) p i„/0 A c/iyfi + ivd/3 A 

Now apply lemma HI □ 
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Lemma 12. If the system is on A m -shell then (J4"5]l and ( |4"6"D /io/d. 
Proof. From 



dA m . , . dA 
IvCb A — h i v d( B A 



d( B ^ d(d( B ) 

/ <9A m \ dA m 
= (-l) PB i v ( B A di^— — —J + i v d( B A 



. dA m \ ,. „ 5A m . „ 5A m 
= -« KCb A + dz„CB A + z„dCfl A 

V o{dQ B )J o(d( B ) o(d( B ) 

/ , <9A m \ „ _ <9A m 

= -a I z„Cb A + £„Cb A 

V d{d( B )J d{d( B ) 

Hence (@S]) follows from (pES}. Likewise from ([27]), ([22]) and d 2 = 

i«CBAd( ) + i v d( B A — h (-1)^ +1 z^Cb A df— — 

= (-l) Pfl ^Cs A d ( — + i v d( B A — 
V dQ? / oQ B 

= 

hence (USD follows from (jSD □ 
Lemma 13. 

£ (/u) a = / C v a + df A for / G TA°M , v E TTM and a e TA P M (118) 

Proof. From Cartan's identity 

£(fv) a = i(f v )da + di(f v )da> = fi v da + d(fi v a) 

= fi v da + df A z^a + fdi v a = f C v a + df A i v a 

□ 

Lemma 14. (15 9p and (|60|) are equivalent if Q e rA™M zs non-vanishing. 

Furthermore if (1601) zs inte /or one non-vanishing Q it is true for all non-vanishing 
Q e TA n M. 

Proof. 

T(u, v)Q = g(T(u),v)Q = i-r^vQ = v A ir(u)^ 



hence (1591) and (1601) are equivalent if f2 6 rA n M is non- vanishing. 
Substituting ft -»■ = where J e TA°M is non- vanishing then 

v A ir(u)& — u A ir(v)& = </(w A ir(u)Q — u A iq-^Vt) = 



□ 



29 



Lemma 15. Given a map r : TTM -» TA^-^M, r : v ^ r v , let T : TTM TTM be 

defined with respect to the non vanishing top form Q via r v = ij-Vt. 

Using a coordinate frame (x 1 , . . . , x n ) with T(v) a = dx a (T(v)) and Q = dx 1 A • ■ ■ A dx n 
implies 

dx a A r v = T{v) a Q (119) 
T(v) a = i dn ---i dl (dx a Ar v ) (120) 

and 

dr v = d a (T(v) a )Q (121) 
Using an orthonormal coframe {e 1 , . . . , e n } with l~(v) a = e a (T(v)) and Q = *1 then 

e a A r v = T(v) a n (122) 

and 

T{v) a = *~\dx a At„) (123) 
Proof. Using the coordinate frame r v = ir(v)^ = T{v) a iQ a VL 

dx a A t v = T{v) b dx a A i db Q = T{v) b = T{v) a Q 

since clearly 



dx a A i db Vt = dx a A z^cfe 1 A • • ■ A dx n 

= {-lf^dx A dx 1 A • ■ ■ A dx*' 1 A dx b+l A • ■ ■ A dx n 



= dx 1 A • • ■ A dx*" 1 A dx a A dx b+1 A ■ ■ • A dx n = 5%dx l A ■ ■ • A dx n 
Hence f lll9|) and since z a „ • • • i^cfe 1 A • • • A ab n = 1 (I120p . Also 

dr v = d(i nv) n) = d(T(v) a i da n) = d(T(v) a ) a i da n + T(v) a d(i da n) 
= d b (T(v) a ) dx b a i d n = d a (T(v) a )n 

Equation f II 2 2 j) is proved similar to f)119p and f)123p is trivial. □ 
Lemma 16. (EE) holds 



Proof. Using a g-orthonormal frame X a and dual coframe e a , set rj ab = g(X a ,X b ) £ 
and T H = T ab e a ® e b . Then from 



AA m 

T ab * 1 = *T w (X a , X fe ) = 2 — — i(e a <g> e b ) 

Ag 



hence 

AA 



. = T ab (*l) ® X a ® X b 
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From the algebraic symmetry of T n , T ab = T ba . 
From ([HID wit h v = v a X a 

= *(T>, -)) = *((T ab e a ® e b )(v, -)) = v a T ab * e b 

thus using ([3"9"[) and f[3"Tj) with = *1 yields 

/AA m \ / \ 

2 A , 9) = A (T ab (*l) <g> X a ® X 6 , e c ® e c J 

= A ((X 6 : e c )(*l) ® (T a6 X a ), e c ) + A ((T ab X a : e c )(*l) ® X 6 , e c ) 

= A ((*1) ® (T ac X a ), e c ) + A (r/a C T afe (^l) ® X 6 , e c ) 

= e c (v)w * 1 + e c (^)T a c ^ Xa * 1 = v c T ac * e a + ^ c T a c * e a 
= 2v c T ac * e a = 2t v 

□ 

Lemma 17. (EI]) /ioWs 

Proof. Again using an orthonormal frame and setting T H = T a i,e a ® e 6 and u = i? a X a as 
in lemma [16] and using ([39]) , ( [37]) , the algebraic symmetry of T n and ([90]) yields 

/ AA m \ / \ 
2B„ , g) = B v (T ab (*l) ® X a <g> X 6 , e c ® e c J 

= S„((X 6 : e c )(*l) ® (T afe X a ), e c ) + B„ ((T af X a : e c )(*l) ® X 6 , e c ) 
= B„((*l) (8) (T ac X a ),e c ) + i3„(V c 6 (*l) ®X 6 ,e c ) 

= (_ 1 )n-l T ao^ ea) A -^ ec _ d[T ac ^ ^ A 

+ (-l)™- 1 ^ 6 ^) A z„rfe c - d{T b * e 6 ) A i v e c 
= 2(i v de c A T ac (*e a ) - d(T ac * e a ) A ^e c ) 

= -2(v c d(r£ ) - i v de c A r« ) = -2(Dt h ) v 

□ 

The remainder of this appendix is to derive ( [52]) in the text. Since it is convenient 
not to use the summation convention for multiindex summations, the following lemmas 
require the introduction of some new notation. 

Given positive integers n and r, where r < n then denote Inc(r, n) as the set of lists 
of r elements, which are increasing 

Inc(r,n) - {/ : {1, . . . , r} ->■ {1, . . . , n} | h < I 2 < ■ ■ ■ < I r } 

We use the notation / 6 Inc(r, n) with len(J) = r. For example Inc (2, 3) = {[1, 2], [1, 3], [2, 3]}. 
The sum over the set Inc(r, n) is then written 

n n n 

E = E =E E ••• E 

7Glnc(r,n) l<h<h<-<I r <n Ii=l I 2 =h+1 I r =I r -i+l 
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Given positive integers n and r, Free(r, n) is the set of lists of r elements not necessarily 
increasing 

Free(r,n) = {i:{l,...,r}->{l,...,ra}} 

We use the notation / G Free(r, n). For example 

Pree(2, 3) = {[1,1], [1,2], [1,3], [2,1], [2, 2], [2, 3], [3,1], [3, 2], [3, 3]} 
The sum over the set Free(r, n) is then written 

n n n 

E = E =EE-E 

/GPree(r,n) l</i,...,/ r <n 7 1= 1 7 2 =1 / r =l 

Let {e a , a = 1, . . . , n} be a coframe (not necessarily orthonormal) and X a its dual. Let 

e 1 = e h A • • • A e Ir 

and let z a = %x a and 

ii = ii r ---ih (124) 

Clearly for antisymmetric then 

E ^ e/ = ^[ £ ^ e7 ( 125 ) 

ielnc(r,n) /eFree(r,n) 

For free list / G Free(r, n) let e(J) is the signature of / and Inc(J) is the increasing 
form of /. Here e(J) = if / contains repeated indices. Thus 

e 1 ' = e(/)e Inc(/) (126) 

For free lists /, J G Free(r, n) let 

/ _ f e(I)e(J) if Inc(l) = Inc(J) 
J 1 Otherwise 

Thus for increasing lists /, J G Inc(r, n) then 5j = 1 if / = J and 55 = if / ^ J. 

We use concatenation to represent the combining of lists, so that if / G Free(r, n) and 
J G Free(s, n) then I J G Free(r + s, n) in the natural way, i.e. 



Ifi if A* < r 

i M -r if A* > r 



Likewise if a G {1, . . . , n} and / G Inc(r, n) then a/ G Free(r + 1, n). We use the backslash 
to represent the removal of an element from a list. That is for J G Free(r, n) and 1 < s < r 
then J\J S G Free(r — 1, n) is given by 



(M)m 



J M if /i < s 

if A* > s 
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Lemma 18. 

e a Ai a a=pa for a = TA P M (127) 

Proof. Clearly true for p = 0. By induction, assume true for a G TA P M let f3 G TA 1 M. 

e a A i a (/3 A a) = e a A (i a /3 A a) — e a A (/? A 2 a a) = /3 A a + (3 A e a A z a a 
= (p+ 1)P A a 



□ 



Lemma 19. 



7elnc(r,n) JeFree(r,n) 

Proof. Clearly the first two expression are equivalent. 
Consider r = 1 then 

i a (e a A a) = 5^a — e a A i a a = (n — p)a 
Then the left hand side of (I128p becomes 



' V ^ \it ■ ■ ■ if (e fl A ■ ■ • A e ir A a) 



J2 hi** 1 A «) 

j£ Free (r,n) 



-j- (n - (p + r - 1))^ • • • if 2 (e^ 2 A • • ■ A e' r A a) 

■ ■ ■ — —(n — (p + r — 1)) • ■ • (n — p)a 

1 / ,v / \ (n — p)! 

— (n — p — r + 1) ■ ■ • (n — p)a = — -a 

r! r\[n — p — r)\ 



□ 



Corollary 20. 7/deg(e / A a) = n then 



ii(e I Aa)=a (129) 

/£lnc(r,n) 

Proof. Since n — p = r. □ 
Lemma 21. For increasing lists I, J G Inc(r, n) and f2 G TA n M then 

e 1 Ai J n = nS I J (130) 
Proof. In the proof of lemma [TS] we show 

e a A i b Q = 

So 

e J A ijVt = e h A ■ ■ ■ A e Ir A i Jr ■ ■ ■ i^Q = e h A ■ ■ ■ A e^" 1 A i Jr _ Y ■ ■ ■ ij^i 



ns^s'r) ■■■5 I j i = n<rt 



□ 
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Corollary 22. For free lists I, J G Free(r, n) and Vt G rA n M i/ien 

e 1 A ijVl = QSj (131) 

Proof. 

e 1 A i,n = 6(/)e(J)e Inc ( f ) A z Inc(i) Q = e(J)e(J)a£$n = <5jfi 



□ 



Lemma 23. Let I G Inc(r, ?i) and J G Inc(r + l,n) and Q G TA n M t/jen 

r+l 

e ; A^ = ^(-l) s -% Js ^ (132) 

s=l 

Proof. First observe that for the left hand side of (11321) to be non zero then I C J. Since 
J and J are increasing then J r+ \ ^ J s for s = 1, . . . , r — 1. Thus there are two cases, 
either J r+ i = I r or J r+ i does not equal any I s . I.e. 



ij r 



/ = (-l) r -%: +i e IVr (133) 



By induction on r. Clearly true when / G lnc(0, n). Assume true when / G Inc(r — 1, n) 
then 

e 1 A ijfi = e 1 A ij r+1 ij\ Jr+1 £l 

= (-irz Jr+1 (e 7 A Oy^fl) + (-l)'"^^ A ijy^n 

from (JI30D and ([133]) 

= (-i)^W^ + a i JVr+1 n 

from induction hypothesis 

r 

= (-ir<W^ + ( B-^^C+a^^) 

r 

= (-±YSj\j r+1 u r+ ^ + E( \r '^,o..<> 

s=l 

r+l 

□ 

Lemma 24. Le£ V J be a multiindex object where J G Inc(r + l,n) and /e£ / G Inc(r, n) 
and tt G TA"M then 

r+l rt 

E V ' X>l)*~%J.V.n = E V Inc(a/) e(a/)z a fi (134) 

J€lnc(r+l,n) s=l a=l 
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Proof. The only values of J such that = 1 are when J = Inc(aJ) for some a. In this 



case one value of s is 5js Js = 1 and 



Hence result 



\s-l 



e(al) 



□ 



Corollary 25. Let V J be an antisymmetric multiindex object where J G Free(r + l,n) 
and let I G Inc(r, n) and f2 G TA ra M i/ien 



r+1 



JGlnc(r+l,ra) *=1 



(135) 



a=l 



Proof. Follows since 



V Inc(a/) e(aJ) = V a/ 



□ 



Corollary 26. Lei V J and a/ be antisymmetric multiindex objects where J G Free(r+1, n) 
and I G Free(r, n) and let VL G TA n M then 



/Glnc(r,n) JElnc(r+l,Ci) /Glnc(r,n) o=l 



(136) 



Lemma 27. ([52]) /io/ds. 
Proof. From (J49]) 



Let deg(dCe) = r + 1 and d£s — J2ieinc(r+i,n)(d(B)idx I and J G Inc(r + l,n) then from 
(ESD 



da; J A 



dA m d 



d(d( B ) de 
d 
de 



£=0 



A m (...,d( B + edx u 



OS 



s=o o{d( B )j 



n 



Hence from (11291) 



dS 



Jelnc(r+l,n) v v 



— ijtt= V ij(dx J A 

B)J T ^r,i s V 



dA m \ dA 1 



Jelnc(r+l,n) 
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From (1471) 



M da = ia A m -J2£ da (B 



A 



aA 1 



B=l 



d(d( B ) 



A 
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B=l 7glnc(r,n) Jelnc(r+l,n) 
Q 

«si n-E E E 5 «(c 

-8=1 Ielnc(r,n) Jelnc(r+l,n) 

Q 



9« 



bJj 



55 ,7 • ^ 



<9(dC 



b; j 



<ts^-E E d ^i) 

B=l /Glnc(r,n) 
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Hence 



Thus 



giving 



Now since 



then 



and 



I.e. (E2D- 



from 



BiW 



= 



\ B=l /elnc(r,n) 

= (rV-^ + E E ^ 

^ B=l /eInc(V,n) 
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d(d( 



Q 



rV = <&s-E E fl -(CBz) 

B=l 7glnc(r,n) 



B 
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9(dC 



^Cb= E ^(Cb/^)= E ^Bidx bI = E e(6J)9 fe CB/^ Inc(W) 

/Glnc(r,n) /Glnc(r,n) 7glnc(r,n) 



95 



95 



e(blf 



95 



95 



d(d( B ) bI d(d( B )i nc{bI) d(d b ( BI ) d(d b ( BI 
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r M a b = s b a s-J2 E d ^i) 
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B=l /Glnc(r,n) 
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B Computational method for calculating Gateaux deriva- 
tives 

In any local frame {X a } with dual co- frame {e b } let A m (Z, . . .) e TA n M depend on the 
tensor Z = Z a b X a ®e b . Relative to any chosen top-form f2 write A m (Z, . . .) = S(Zg, . . .)Q 
where S G TA°M and is independent of Z. Thus S(Z a b , . . .) depends on the n 2 variables 
Z a b with n 2 derivatives -§§zr- In this case 



AAm :y_ d 
AZ ' ~ de 



d 



e=o cfe 



95 



A m (z + eY,...) = - s{z\ + £ y a 6 , ...)n = ^—Y a b n 



e=0 y °' ' az- 



6 



(95 

-n <g> e a <g> x 6 ;f 



hence 

AA m dS 



il(g)e a 0X b 



AZ dZ a h 

For example if S is such that A m (Z, a, (3, u, v) = Z(a, v)Z(/3, u)Q for any 1-form fields 
a, (3 e TA 1 M and vector fields u,v E TTM then 

= Z a b Z c d a a v b p c u d 

and it follows from the definition that 
AA m 

" ,i"7C O ,0\ r, a V i „,d/3 „M r 7 c 



a a v a p c u b Z c d Q ®e a ®X b + a c v d (3 a u b Z c d Q ® e a ® X b 
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